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Abstract. A theoretical investigation of the bound-free electron-positron pair production in relativistic 
heavy ion collisions is presented. Special attention is paid to the positrons emitted under large angles with 
respect to the beam direction. The measurement of these positrons in coincidence with the down-charged 
ions is in principle feasible by LHC experiments. In order to provide reliable estimates for such measure- 
ments, we employ the equivalent photon approximation together with the Sauter approach and derive 
simple analytic expressions for the differential pair-production cross section, which compare favorably to 
the results of available numerical calculations. Based on the analytic expressions, detailed calculations are 
performed for collisions of bare Pb 82+ ions, taking typical experimental conditions of the LHC experiments 
into account. We find that the expected count rate strongly depends on the experimental parameters and 
may be significantly enhanced by increasing the positron-detector acceptance cone. 
Key words. Bound-free pair production, Relativistic heavy-ion collisions, Virtual photons 
PACS 25.75.Dw, 12.20.Ds, 25.75.-q, 03.65.Pm 

PACS. XX.XX.XX No PACS code given 



1 Introduction 

Owing to recent advances in heavy-ion accelerators, an 
increasing interest arises in exploring electromagnetic pro- 
cesses accompanying ion collisions. One of the dominant 
' processes is electron-positron pair production; pertinent 
5h '. cross sections are large. For ultra-relativistic collisions of 
C3 • two bare lead ions (Pb 82+ ), for example, the cross section 
of the creation of a free e + e~ pair may reach hundreds of 
kilobarns [1-5]. Both free-free as well as bound-free pair 
production can take place in which the electron is captured 
by one of the projectiles resulting, thus, in the formation 
of a hydrogen-like ion, 



Z 1 + Z 2 -> Z x + e+ + {Z 2 + e") . 



(1) 



where the bound system is denoted by round brackets. 
Even though this (bound-free) process is usually orders 
of magnitude less probable than free-free pair production, 
its investigation is of great importance not only for a bet- 
ter understanding of the physics of extraordinary strong 
electromagnetic fields but also for the development and 
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operation of novel collider facilities. In free-free e + e~ pro- 
duction, the scattered nuclei lose only a very small fraction 
of their energy and acquire tiny scattering angles, and thus 
do not leave the beam. In contrast, if one of the colliding 
ions captures an electron [cf. Eq. (1)], it changes its charge 
state and is bent out from the beam. The corresponding 
cross section is rather large, about 100 barn, and the re- 
action (1) is one of the important processes which limit 
the luminosity of colliders. Besides, the secondary beams 
of down-charged ions emerging from the collision point hit 
beam-pipe and deposit a considerable portion of energy at 
a small spot, which may in turn lead to the quenching of 
superconducting magnets [6-9]. 

Because of its fundamental and practical importance, 
the bound-free pair production has been in the focus of 
intense research over the past years. A series of exper- 
iments have been performed, for example, at the CERN 
Super Proton Synchrotron (SPS) to analyze the total cross 
sections of the process for ultra-relativistic collisions of 
highly-charged Pb ions with solid-state and gas targets 
[10, 11]. These experimental findings are currently under- 
stood based on theoretical predictions of relativistic Dirac 
theory [12, 13]. In contrast to the total rates, much less at- 
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tention has been paid until now to the differential bound- 
free pair production cross sections. However, the angle- 
resolved analysis of the positron emission is of definite 
interest since it allows to probe a parameter range (of the 
e + e~ process) which is otherwise not accessible in total 
cross section studies. Namely, while the total probabilities 
arise from a region where the transverse momentum of 
the produced positron is small, p+±_ < m e , large values of 
p+±_ may significantly contribute to the angle-differential 
rates. 



An understanding of the differential pair production 
properties is also significant for the analysis of future 
experiments at the LHC facility. In these experiments, 
the down-charged ions can be measured in coincidence 
with the emitted positrons in order to reduce background 
events and clearly distinguish the process (1). However, 
since the central detector is likely to be placed at a rather 
large angle with respect to the beam direction it will 
only detect those positrons, whose transverse momenta 
are much greater than the electron mass, p + j_ 3> m e . An 
estimate of the yield of these positrons is highly desirable 
for estimating the feasibility of future measurements; the 
task has not been fully tackled until now to the best of 
our knowledge. 

In this contribution, therefore, we present a theoretical 
study of bound-free pair production with a special empha- 
sis on the differential cross sections. The exact computa- 
tions of these cross sections for the ultra-relativistic ion 
collisions and large positron momenta are very laborious, 
and we shall develop an approximate method that allows 
for an analytic treatment. Before discussing the deriva- 
tion of the analytic expressions, we introduce in Section 2 
the notation used in the manuscript and estimate typical 
values of kinematic parameters corresponding to the LHC 
experiments. By using these estimates and the equivalent 
photon approximation (EPA) we express in Sect. 3.1 the 
differential probabilities of the process (1) in terms of the 
pair p/io£o-production cross sections. The evaluation of 
these cross sections within the framework of the Sauter 
approximation is discussed in Sects. 3.2 and 3.3. The va- 
lidity of our EPA-Sauter model is tested numerically in 
Sects. 4.1 and 4.2 where we calculate the total as well as 
differential pair production rates and compare them with 
the predictions of an exact relativistic theory. Based on 
the results of this test, we employ the EPA-Sauter ap- 
proach in order to estimate the positron yield for ultra- 
relativistic collisions of bare Pb 82+ as relevant to the LHC 
studies. In Sect. 4.2 we focus on two experimental scenar- 
ios: In the first scenario, only positrons with very large 
transverse momenta, p+± ~ 1 GeV are "seen" by the de- 
tectors, while much smaller momenta, p+±_ ~ 0.05 GeV, 
are considered in the second case. We find that the reduc- 
tion of the (accepted) minimum momentum p+± leads to 
an increase of the positron count rate from just a single 
event in 67 days to about 10 events per hour. We conclude 
with a brief summary in Sect. 5. Through the paper, we 
use relativistic units with c = h = 1 and a ~ 1/137. 



2 Notations and kinematic parameters 

Let us first recall basic notations and assumptions used 
throughout this paper. The initial state of the overall sys- 
tem is given by two bare ions of charges Z\ = Z2 = 
Z and masses M\=Mi=M, moving towards each other 
with 4— momenta = (£1,27 -Pi, 2) and corresponding 
Lorentz-factors 71 = 72 = 7. Being defined in the col- 
lider frame, these kinematic parameters are convenient for 
the description of experimental results. However, the the- 
oretical study of atomic processes accompanying ion-ion 
collisions is most conveniently done in the rest frame of 
that particular ion which finally becomes hydrogen-like. 
For definiteness, we shall assume electron to be captured 
by the "second" nucleus in whose rest frame the "first" 
nucleus moves with a Lorentz factor 7^ = 27 2 — 1. By 
adopting 7 = 1500, which is the typical value for LHC 
experiments on Pb-Pb (Z\ = Z2 = 82) collisions, we find 
7i = 4.5 -10 6 . This value is used in all numerical estimates 
below. 

Bound-free pair production in energetic ion collisions 
can be uniquely detected experimentally by measuring 
the down-charged ion in coincidence with the emitted 
positron. In the set-up of the typical LHC experimen- 
tal arrangement, the positron detector is placed at rather 
large angles with respect to the incident beam direction. 
Therefore, only positrons with a transverse momentum 



P+± > Pram > m 



(2) 



can be observed in a particular experiment. In Eq. (2), 
m = m e is the electron mass and momenta are given in the 
collider frame. In the same frame, we define the positron's 
rapidity: 



1 , £++ P+z 

v+ = o ln 

2 e+ - P+z 



-ln[tan(±0+)] , (3) 



where e+ and p+ z denote the energy and longitudinal 
(along the beam axis) momentum of the positron and 8 + 
its emission angle. These quantities are directly related to 
the transverse momentum, as follows, 



P+J 



P+z 



Pa 



tan 9 j 



(4) 



Equation (2) immediately implies restrictions both on the 
rapidity, 

- Vmin < V+ < 2/min , (5) 

as well as on the positron angle, 

n , 6 min = 2 arctan (e" Wmin ) . (6) 



9min < d+ < Tr- 



hi bound-free pair production, certain limitations are 
imposed on the electron's kinematic parameters. As the 
electron is produced in bound ionic states, its relativistic 
Lorentz factor and scattering angle match those of the 
"second" nucleus. Indeed, for ultra-relativistic collisions 
of heavy high-Z ions, the corresponding scattering angle 
should be small, 



< 



E-2 



P+±_ 

jM 



10" 



(7) 
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Fig. 1. Schematic Feynman diagram for the bound-free pair 
production in a heavy-nucleus collision. Heavy colliding nuclei 
are represented by thick lines while thinner lines correspond 
to light fermions (electron and positron) . The trajectory of the 
produced electron is aligned with that of the second nucleus. 



Here, all parameters are given in the collider frame. It fol- 
lows from Eq. (7) and the previous discussion, that the 
electron's momentum is almost parallel to the initial di- 
rection of propagation of the "second" nucleus, chosen as 
the z-axis: 



M«0, 



V- 



OS) 



Using these expressions and Eqs. (2)-(4) for positron 
emission, we can find the four-momentum p = (e,p) = 
(e+ + £-, p+ + p_) and the invariant mass 



W = 



> 



V 2 P+P- = V 2 P+-L m l tan {\ ^ 



2pminTO7 tan 



(9) 



of the final e + e~ system. 

Our theoretical analysis of the bound-free pair pro- 
duction is based on a quantum electrodynamic (QED) ap- 
proach. In lowest order, the process (1) is described by the 
diagram in Fig. 1. The exchange of two virtual photons, 
with four momenta q\ y i = (wi,2, q.1,2)) is required to pro- 
duce the e + e~ pair. For the present study, it important to 
know the virtuality Q 2 of these photons which measures 
how far they are off the mass shells. For the first gamma 
quantum with a 4-momentum q\ = Pi — P{, for example, 
we write 



where we use the fact that 



M 2 

+ E', - 

2E 1 1 



M 2 



2E' 



The virtual photons are "emitted" 
counter-propagating nuclei, and we finally obtain 



<l2z 



-L02 ■ 



(10) 

(11a) 
(lib) 
by 

(12) 



Together with the four-momentum conservation relations 
uji + u>2 — £ and q\ z + q2 Z = Pz and by employing Eqs. (4) 
and (8), this allows us to derive the energy of the first 
virtual photon, 



1 



(e-Pz) = o ( £ + -P+z + 



E- 



P-z) 



sin ( 



tan( 



tan ( i 9 



(13) 



as a function of the transverse momentum and emission 
angle of positron. Employing the restrictions given in 
Eqs. (2) and (6) on these quantities, the minimum allowed 
value of bj\ can be found, 



■ r 1 1 

mmjcji} = -pmm tan 



(14) 



Before we proceed to the theoretical treatment of the 
pair production, let us first compute some typical 
values for the basic kinematic parameters from above. As- 
sumptions about the magnitude of these values will be 
required later in order to perform a number of approxi- 
mations that greatly simplify the calculations of the cross 
sections. As we have seen already, the ranges of positron 
as well as virtual photon parameters are crucially depend 
on the particular experimental set-up. Here, two distinct 
experimental "scenarios" are considered. 

In the first scenario, which relies on a typical LHC- 
detector acceptance, the transverse momentum and ra- 
pidity of the positrons are restricted by the conditions: 



Pn 



1 GeV. 



Vu 



1 



For these parameters, we find 

0mi„ = 40° , 

W > 0.75 GeV, 
min{wi} = 0.18 GeV, 



(15) 

(16a) 
(16b) 
(16c) 



for the minimum values of the positron emission angle 
# m i n , the virtual photon energy min{wi} as well as for the 
e + e~ invariant mass W , respectively. 

In the second scenario, we assume smaller transverse 
momenta of the positron and a greater rapidity: 

Pram = 0.05 GeV , y min = 1.5, (17) 

which results in significantly different parameters: 



$min — 25 , 

W > 0.13 GeV, 
min{wi} = 5.6 McV. 



(18a) 
(18b) 
(18c) 



As we see below in Section 4.2, the prospects for a unique 
detection of bound-free pair production seem much more 
favorable for the parameters listed in Eq. (18) than for the 
parameters in scenario (16). 
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3 Theoretical background 

3.1 Equivalent photon approximation 

Since its introduction by Fermi [14] and further devel- 
opment by von Wcizsacker and Williams, the equivalent 
photon approximation (EPA) has been successfully ap- 
plied to the description of a large number of electromag- 
netic processes induced in collisions of charged particles. 
Usually, one physically justifies this method based on the 
observation that the electromagnetic field of a fast mov- 
ing charge becomes almost transverse, and the electric and 
magnetic fields both have about equal strengths [15]. In 
the observer (laboratory) frame, therefore, the projectile's 
fields can be seen as those of the pulse of a plane, linearly 
polarized wave. The frequency spectrum of such a pulse is 
calculated within a semiclassical model and used in order 
to evaluate collision-induced cross sections [16]. 

Even though the conventional EPA is found very use- 
ful for the theoretical analysis of ion collisions, here we 
will employ an alternative approach which — in the ultra- 
relativistic domain — appears to be more transparent and 
rigorous and allows for a simple estimate of its accuracy. 
Our approach exploits the EPA as an approximate method 
for calculating a Feynman diagram for the corresponding 
process and uses the fact that the virtual photons in the 
diagram are close to the mass shell — for details, see the re- 
view [1]. The diagram for bound-free pair production (1) 
is displayed in Fig. 1, where the two thick lines represent 
the colliding nuclei and the double-line arrow just refers 
to the residual hydrogen-like ion. Inspecting this diagram, 
we may interpret the bound-free pair production as being 
due to the interaction of the second (lower) nucleus with 
the virtual (or equivalent) photon with the energy uji and 
virtuality Q 2 = —q\ emitted by the first nucleus. Thus, 
the theoretical analysis of the e + e~ creation in energetic 
heavy-ion collisions can be traced back to the virtual pro- 
cess: 

1 * + Z 2 ^e+ + (Z 2 + e-) ls , (19) 

where the electron is captured into the ground ionic state. 
The differential cross section dazz of the process (1) can 



be expressed, therefore, in terms of the cross sections aE 



and cr^z f° r ^ ne (virtual) pair production (19) as: 



dazz = dux da. 



7*Z 



dna da. 



7*Z 



(20) 



In this expression, we denote by driT and dns the numbers 
of the equivalent transverse T and scalar (or longitudinal) 
S photons and, moreover, neglect the interaction between 
the emitted positron and the first nucleus. 

As seen from Eq. (20), any further evaluation of 
dazz requires the knowledge of both of the pair photo- 

T S 

production cross sections da^ z and infinitesimal numbers 
driT.s of equivalent photons. Below we will evaluate these 
quantities in the rest frame of the second ion, which finally 
constitutes a hydrogen-like ion. The energy of the equiv- 
alent photon, emitted by the first nucleus (cf. Fig. 1), is 
given in such a frame as: 



qiP 2 
M 



27W1 



For the case of an LHC experiment, the minimum value of 
this energy, min{wi} = 27min{tJi}, is huge: it amounts 
to about 550 GeV and 17 GeV for the first (16) and second 
(18) scenarios, respectively. 

In ultra-relativistic nuclear collisions, the main contri- 
bution to the pair production cross section (20) is given 
by a wide range of the photon virtualities (10) ranging 
from the small minimum value 



I 2 . 

J mm 



7i 



(22) 



to some Q m ax which is derived below from the analysis 

T S 

of driT.s an d the differential cross sections da^ z . In such 
a region, the numbers of equivalent photons are given by 
(see Appendix D of Ref. [1]): 



dn s (u> L , Q 2 ) ~ dn T (u L , Q 2 ) 
_ Zfa du L dQ 2 ( \ Q 



TT LO L Q 2 



1 



2 

min 



Q 2 



F 2 (Q 2 



(23) 



where F(Q 2 ) is the form factor of the first nucleus. 

Integrating Eq. (23) over the interval [Q min , Q max ]i we 
find the overall numbers dns(wi) and dnT^h) of the 
scalar and transverse virtual photons which contribute to 
the production process (1). For the g?77/t(wl) we obtain: 



dn T (wL) 



Z 2 a du>L 

TT UJ L 

" »« dCf 

Q 2 



) 2 . 

' min 



Q 2 



F 2 {Q 2 ). (24) 



The evaluation of this integral is discussed in detail in 
Ref. [17] and employs the Q 2 -behavior of the (square of) 
the nuclear form factor. For example, the function F 2 (Q 2 ) 
drops quickly when the virtuality of the photon becomes 
greater than the squared inverse electromagnetic radius of 
the nucleus, Q 2 > 1/R 2 . Therefore, if C max » 1/-R 2 , one 
can simply extend the upper limit of the integration in 
Eq. (24) to infinity and obtain: 



where 



dn T (uj L ) 



Z 2 a diVL 



g(u> L R/jl) 



°°dy (. x 2 



^ (1--) F 2 (y/R 2 ) 

v V v 



(25) 



(26) 



In particular, it was found that for small values of x = 
wlR/jl,i which correspond to the range of parameters 
considered in the present article, g{x) can be approxi- 
mated with an accuracy better than 1% as 



SO) = In ^2 



Cq ; 



(27) 



with C = 0.163 for Pb and C = 0.166 for Au. 

Equations (25) and (27) provide an approximation to 
(21) the photon number dnr(wi) in the regime where Q max S> 
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1/ R 2 . A simple estimate of the driT^L) can be obtained 
also for the Q^ ax <C 1/R 2 . In the latter case, the nucleus 
can be treated as pointlike with the F(Q 2 ) = 1. By in- 
serting this form factor into Eq. (24) we find: 



and do not exhibit logarithmic enhancement. Therefore, 
the accuracy of the EPA as determined by the relative 
order of the neglected terms is determined by the factor 



1 



dn T (wi) 



Z\ a duiL 



UJ L 



In 



ilQ 



2 

max 



m 



(28) 



{RW) 2 i 



< 1%. 



(34) 



Having clarified the behavior of the photon numbers 
driT.s, we now turn to the question of how the cross sec- 
tions dt7^f z (ujL, Q 2 ), which also enter Eq. (20), depend 
on the virtuality Q 2 . In order to perform this analysis, we 
consider pair production in the fields of virtual and real 
photons. This 7*7 — > e + e~ process is known [18] to ex- 
hibit very similar behavior as (19), but can be described 
by simple analytic expressions given in Appendix E of 
Ref. [1]. Namely, while the pair production cross sections 
a^(W 2 ,Q 2 ) drop drastically at Q 2 > W 2 , the following 
estimates are valid for smaller values Q 2 <C W 2 , 



For the total cross section, the main region corresponds 
to W 2 ~ 4m 2 <C 1/R 2 , and we should use the expres- 
sion (28) for the number of equivalent photons dnT(wi). It 
contains a large Weizsackcr-Williams logarithm ~ Iwy^. 
On the other hand, the omitted items are of the order of 
unity, 

rim2 dQ 2 Q2_ 
Q 2 W 2 



1 



(35) 



and again have no logarithmic enhancement. Therefore, 
the accuracy of EPA in this case is of the order of 



1 



7*7 - T/J/2 CT 77 ' 



Hil) 



3%. 



(36) 



7*7 



y 77 



l + O 



w 2 



(29) 



3.2 Bound-free pair photo-production in the Sauter 
approximation 



where 



11 (W ) describes the e" 
by two real photons. 

We are now ready to sum up results of the above anal- 
ysis and to evaluate the cross sections for the bound-free 
pair production. For the photon virtualities in the inter- 
val [Q^in, Qmax] ( see discussion of Eq. (22)) we can neglect 
the contribution of the scalar photon and approximate the 
cross section of the pair production by the transverse pho- 
ton da^ z by its value dcr^ z on the mass shell (i.e. when 
Q 2 = 0). Under these assumptions, Eq. (20) simplifies to 



pair production Finally, as seen from Eq. (30), the computation of 



da 



EPA 
ZZ 



dn T {w L ) da lZ (<jj l , p + ± ) 



(30) 



The accuracy of such an approximation should be verified 
for two kinematic regimes. If W 2 ^> 1/R 2 , which corre- 
sponds to the experimental conditions of an LHC detector, 
the photon's virtuality is restricted by the nuclear form 
factor to 

Q 2 < ^ « W 2 (31) 

and we should use the expressions (24)-(27) for the num- 
ber of the equivalent photons dnxi^L)- iFrom Eq. (27), 
we obtain the large Weizsacker-Williams logarithm £ — 
In [^/(Rojl) 2 ] that originates from the integral 



V R2 dQ 2 
. Q 2 



On the other hand, the omitted terms are small, 



(32) 



V* 2 dQ 2 Q 2 1 

Qlnn 



} 2 W 2 R 2 Q 2 



(33) 



dcrf P Z A within the EPA can be traced back to the differ- 
ential cross section of the bound-free pair production fol- 
lowing (real rather than virtual) photon impact on a bare 
nucleus: 

1 + Z + {Z + e-) ls . (37) 

During the past two decades, this process has been stud- 
ied in detail with a special emphasis on high-Z ions. The 
total as well as differential cross sections have been evalu- 
ated within a relativistic framework in Rcfs. [19, 20]. The 
fully relativistic calculations, based on the partial-wave 
representation of the Dirac continuum states, were found 
to provide accurate predictions in the near-threshold re- 
gion but face well-known problems connected with the 
slow convergence of the multipole expansions when the 
photon energy increases. Some approximate methods have 
to be used, therefore, in order to calculate the cross sec- 
tions of the process (37) for very high photon energies 
that correspond to large transfer momenta in ion-ion col- 
lisions of the type (1). In the present work, the calcula- 
tions are based on the Sauter approximation (SA) which 
was originally derived for the atomic photoeffect by as- 
suming ultra-relativistic electron energies and disregard- 
ing terms of relative order a 2 Z 2 in an aZ-expansion of 
the transition amplitude [21-23]. The general formulas of 
such an approximation, when applied to the bound-free 
pair photo-production, are obtained below, while their va- 
lidity for the high-Z is discussed later in Section 4.2. 

The differential cross section for the process (37) can 
easily be deduced from Sautcr's relativistic formula for 
the photoeffect or photorecombination, in the notation of 
Eq. (37) as given, for example, in Eq. (57.8) of Ref. [24] 
upon replacing the electron four-momentum p in the 
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photocffcct with the momentum (— p+) of the outgo- 
ing positron in the bound-free photo-production (crossing 
symmetry). More specifically, this implies the following 
changes: 



7 7l 



v 2 ^ 



ll 



= -v+ , 

1L + 1 
1L 



1 — v cos 6^1 — v + cos , 



(38) 



where v + and are the velocity and the scattering angle 
of the positron in the nuclear rest frame, and u) L = (j L + 
1) m is the photon energy. Performing these substitutions 
and averaging over the incident photon polarization, we 
find 



dQ + 



v + sin ■&+ 



m (7L + 1) (1-v+costf+f 
v\ ( 7L + 2) (1 - v+ costf+) - 2^-^- 



(39) 



where, moreover, an additional factor two was included to 
account for the the different statistical weights of leptons 
in the initial and final states of the photo-effect and the 
pair photo-production. 

Equation (39), obtained within the Sauter approxima- 
tion, describes the emission pattern of positrons created 
in the process (37). Integrating this expression over the 
positron angles $ + , we find the total cross section of pair 
photo production, 



SA 



5^,6 



= 4-7T ■ 



Z°a 



G{il) , 



(40) 



where the function G(7l) is given by: 



(x + iy 

2 



r + 3 



x 

V X 



— In (x + Vx 2 - l) 

': 2 — 1 V ' 



(41) 



One may note that this expression coincides with Eq. (52) 
from Ref. [19] upon the replacement E + = —jl- In the 
high-energy limit, ■yi ^> 1, Eq. (40) reads 



,T SA 

<J 1 Z 



An 



Z 5 a 6 
m 2 "fL 



(42) 



where the dominant contribution is from the first term in 
square brackets in Eq. (41). 

3.3 Evaluation of total and differential cross sections 

In the two previous sections, we have shown how the 
cross section for the process (19) can be expressed in terms 




100 



1000 
Lorentz factor y 



Fig. 2. (Color online) The relative difference (55) between 
rigorous relativistic results and predictions of the scaled EPA- 
Sauter approximation (54) for the total cross sections of 
bound-free pair production. Calculations are performed for 
the energetic collisions between two protons (dash-and-dotted 
line) as well as bare gold Au 79+ (dashed line) and lead Pb 82+ 
(solid line) ions. The EPA-Sauter cross section is given by 
Eqs. (50) and (51) derived using a value of = 4m 2 . 



of pair photo-production cross sections and we have eval- 
uated the latter ones within the framework of the Sauter 
approximation. Now we are ready to conclude this anal- 
ysis and to derive the final expressions that characterize 
the positron emission accompanying bound-free pair pro- 
duction in ultra-relativistic collisions between two ions of 
equal charge: 

Z 1 = Z 2 = Z. (43) 

For the high-energy and large-momentum transfer 
regime, where 7l?ti 3> p+i_ = to7li?+ 3> m, the differ- 
ential cross section (39) simplifies to: 



da^(uj L , P+± ) = 16tt 



Z or 771 m dp^ 



P 3 +j 



(44) 



Inserting this expression into Eq. (30), and using the pho- 
ton number (25), we obtain the cross section 



ZZ rn 2 io L y \ lL ) oj l ( p+± ) 3 1 ' 



iZa) 7 m 3 e y + f p+^Re~ y + , 

—4 9 — — — ) d v+ d P+j 



7 (p. 



+-U 



2~ 



which can be used for an estimate of the positron yield in 
LHC experiments. The function g is defined in Eqs. (25) 
and (26). However, since the central detector in such an 
experiment observes the positrons emitted in a wide range 
of angles, one has to integrate Eq. (45) over the rapidity 
y + in the interval (5) and over the momentum p+±_ > p m in 
to get the partial cross section Aa%\ , which is relevant to 
experimental observation: 



A(T ZZ 



16 {Zaf e y - 
3 m 2 7 



L . 



(46) 
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The parameter L is given by 



2 In 



7 



2y„ 



1.44 



-2y„ 



-2y„ 



(47) 



In addition to the partial positron yield (46), we can 
also employ the Sauter approximation (40) in order to 
estimate the total bound-free pair production in ultra- 
relativistic nuclear collisions. We use Eq. (30) together 
with the number of equivalent photons (28) where we 
adopt <5^ ax = 4m 2 , and obtain (upon integration over 
the photon energy): 



a zz 



dn T (u> L ) cr^|(w L ) 



2m 



{a [ln(4 7 l) -!]-&}, (48) 




1000 
Lorentz factor y 



Fig. 3. (Color online) The same as in Fig. 2, but the EPA- 
Sauter cross section is given by Eqs. (50) and (52) employing 
a value of Q^ax = 1.26 m 2 . 



where the parameters a and b can be evaluated analyti- 
cally, 



n 77 

G(x) = ^- = 0.2175 



137 
630 



G(x) In (x + 1) 



57 707-19180 In 2 
88 200 



d.r 



0.5035. 



(49) 



By taking into account that 7^ is defined in terms of the 
nuclear Lorentz factors as 7l = 2 7 2 — 1 we obtain from 
Eq. (48) the well-known scaling law: 



EPA 
'ZZ 



Alnj-B. 



(50) 



where A and B are predicted to be independent of 7 [25- 
27]. The values of these two parameters, as estimated 
within the equivalent photon approximation, depend on 
the particular choice of the maximal photon virtuality 
Qmax- For example, based on our estimate Q 2 nax = 4m 2 
which was recommended in Rcf. [18], we find: 



.4 



16^ a 



3.479 



(Za) 7 



B = 



[6-a(ln4-0.5)] =2.49 



5.72 Z 7 pb, 

(Za) 7 



m 

= 4.09 Z 7 pb, 



(51) 



which are in the picobarn range (pb). By contrast, the 
different value Q 2 nax = 1.26 m 2 adopted in Ref. [28] yields 



A = 3.479 



(Za) 7 



B = 3.49 



(Za) 7 



(52) 



As seen from Eqs. (51) and (52), while the parameter A 
is not sensitive to the change of the parameter B 

varies by about 40 % which is quite remarkable but still 
within the accuracy of the EPA. 



4 Results and discussion 

4.1 Total cross sections 

Although the main goal of the present paper is to pro- 
vide reliable estimates for the positron yield as might be 
measured by the central detector of an LHC experiment, 
we start calculations with the total pair production cross 
sections. Their analysis, performed for different nuclear 
charges and various collision energies, can help to better 
define the range of validity for our EPA-Sauter model. In 
order to start with such an analysis, let us recall that the 
high-energy limit of the Sauter photo-production cross 
section (42) is known to differ from the exact results by 
some factor: 



exact £ ( r 7\ ~SA 
T 1 Z =J( Z ) (J 1 Z 



(53) 



where f(Z) is a decreasing function of the nuclear charge 
[29-31]. Values for f(Z) are obtained in Rcf. [19] based 
on rigorous relativistic calculations. For example, for the 
e + e~ production in which the electron is captured into the 
ground state of an initially bare ion, these calculations pre- 
dict f(Z) = 0.971, 0.222, 0.216 and 0.196 for Z = 1, 79, 82 
and 92, respectively. 

In the following we shall apply the factor f(Z) to im- 
prove also the accuracy of approximations (46) and (50) 
that describe the bound-free pair production in heavy ion 
collisions. For the total cross section of this process, for 
example, we make the replacement: 



EPA 

a zz 



f(Z) 



EPA 
'ZZ 



f(Z) (Alnj-B) 



As seen from Fig. 2, where the relative difference 
f( z ) ■ °f PA 



5 = 



'zz 



'ZZ 



'zz 



(54) 



(55) 



between the scaled EPA-Sauter cross section and exact re- 
sults <7™^ ct from Ref. [32] is displayed, such a replacement 
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leads to a very good agreement between the predictions 
of the two theories (EPA-Sauter and rigorous relativistic 
calculations [32]). In particular, for energetic collisions of 
two gold (dashed line) and led (solid line) ions the rela- 
tive difference 6 does not exceed 1.5% for Lorentz param- 
eters in the range 100 < 7 < 3000. It is worth mention- 
ing that a slightly worse performance of the EPA-Sauter 
approximation can be observed if in place of the photon 
virtuality <5, 2 nax = 4m 2 , used to derive values (51), the 
Qmax = 1-26 m 2 from Rcf. [28] is adopted. For the latter 
case, the relative difference between the exact and EPA 
predictions may reach 6-^12% in the low-energy domain 
(see Fig. 3). 



4.2 Differential cross sections 

Having briefly discussed computations of the total 
bound-free pair production cross sections, we now esti- 
mate the differential pair production probabilities. In or- 
der to assess the reliability of these predictions, let us first 
verify the accuracy of Eq. (44) for der^, and Eq. (45) for 
da z \. These formulates are based on the Sauter approx- 
imation (39) for the e + e~ photo-production. Since this 
approximation is correct only in the leading order of aZ, 
its validity in the high-Z domain has to be examined. To 
this end, in Fig. 4, we display the differential cross sections 
do° x z ct /df2 + of the process (37) for the photon collision 

with bare led ions Pb 82+ and emitted positron energies 
with jl= 5, 10, 15 and 25. Results of rigorous rel- 
ativistic theory, which employs the standard partial-wave 
decomposition of continuum Dirac wavefunctions [20] , are 
compared with scaled predictions of Eq. (39), obtained as 
follows, 



dJ? + 



f(Z) 



dJ?l 



(56) 



where f(Z) is the same factor as employed in Eq. (53). 
One may note that the previously observed, good perfor- 
mance of this scaling for the total rates (the main con- 
tribution to which comes from forward positron emis- 
sion, -d + < 1/7l), does not automatically guarantee ac- 
curate results for the differential cross sections, especially 
in the region of larger emission angles. Fortunately, the 
results presented in Fig. 4 clearly indicate the validity of 
this scaling. Namely, even though the (scaled) Sauter ap- 
proximation can significantly underestimate exact results 
for forward as well as backward positron emission, the 
range of angles, over which the two theories differ signif- 
icantly, shrinks with increasing positron energy. For ex- 
ample, while for ji = 5, the Sauter approximation repro- 
duces the exact results with the accuracy of about 5 % 
only for the angles 50° < t9+ < 60°, this interval increases 
to 2° < tf+ < 120° for 7 L = 25. These findings justify 
the application of the Sauter approximation together with 
the scaling (44) for the analysis of the pair production in 
ultra-relativistic heavy ion collisions at an LHC experi- 
ment in which large positron transverse momenta, p+±, 
will be monitored. 



Besides the direct numerical "proof" , yet another con- 
firmation of the scaling conjecture (56) was recently re- 
ceived. Namely, the asymptotic behavior of the pair- 
production differential cross section in the region m 2 <C 
^ 2mwi was derived for arbitrary Za in Rcf. [36]. 
This asymptotic reads as 

j exact/. , _ \ fasympfy\ -jg^ ^Op 771 771 2 dp+± 

da z {lul,P+±) = J ' l (Z)lbTr - 5 . 

' m z Wi P+± 



asymp (7) — 



2 (277) 



27-2 



r(27 + i) 



|r(7-ir7)| 3 e- 2 " arccos ", (57) 



where i] = Za and 7 = \f \ — rf- . One obtains from this 
expression /^(Z) = 0.297, 0.287 and 0.257 for Z = 79, 
82 and 92. Therefore, our assumptions about scaled factor 
f(Z) are by about 25 % smaller then the corresponding 
factor / as y m P(Z) for the same values of Z. We note that 
such a discrepancy is comparable to the (expected) ex- 
perimental error and the overall uncertainty of our model 
that arises not only from the use of the Sauter and EPA 
approximations but also from neglecting the electron cap- 
ture into excited ionic states, an effect which will be briefly 
addressed later. 

Based on the accuracy analysis of the scaled Sauter ap- 
proximation (44) for da^ z an d (56) for def^ , we are now 
in the position to compute the positron yield. This yield 
is defined by the partial cross section given in Eq. (46), 
which has to be scaled properly, 



2f(Z)Ao% A z 



(58) 



The factor two accounts for the fact that the electron can 
be captured by either of two nuclei. Assuming a maxi- 
mum luminosity C = 10 27 s~ 1 cm" 2 which can be reached 
in Pb-Pb collision experiments at CERN, and using the 
corresponding sets of kinematic parameters {j/ m in,Pmin} 
(see Sect. 2), we finally obtain from Eqs. (46) and (58) 
that (i) only one event per 67 days is likely to be ob- 
served for the high-transverse-momentum scenario (15), 
whereas (ii) about 16 counts per hour are within the ac- 
ceptance range of the detector within the scenario (17). 

One may note that these estimates were obtained 
based on the assumption that the electron is captured into 
the ground ionic states. The e + e~ pair production accom- 
panied by the formation of excited hydrogen-like ions may 
slightly increase these predictions. However, Refs. [19, 33] 
indicate that the effect of the excited state electron re- 
combination shall not exceed 20 -j- 25 %. 



5 Summary and outlook 

In conclusion, we investigate bound-free pair produc- 
tion in ultra-relativistic heavy ion collisions. Special em- 
phasis in our study is placed on the emission of positrons 
with large transverse momenta which can be observed at 
the LHC collider. In order to estimate the yield of such 
positrons, an approximate method, based on the equiv- 
alent photon approximation and the Sauter theory, has 
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Fig. 4. (Color online) The differential cross section (Z 5 7l)~ 1 da 7 z/dJ7+ (in nanobarn/steradian) of the pair production in the 
collision of photon with bare lead Pb 82+ ion accompanied by the electron capture into the ground ionic state. Results of the 
rigorous relativistic calculations (solid line) are compared to the predictions of the scaled Sauter approximation (dashed line). 
Calculations have been performed in the nucleus rest frame and four positron energies rwyL with jl = 5, 10, 15 and 25. 



been laid out. Within the EPA-Sauter approach, we de- 
rive a simple analytical expression for the differential cross 
section of the pair production process (1). Based on this 
expression, calculations have been performed for the col- 
lisions between two bare lead ions Pb 82+ moving towards 
each other with the Lorentz factor 7 = 1500. At this en- 
ergy, typical for the LHC, we consider two scenarios that 
correspond to different detectors set-ups. While for the 
first scenario (15) the number of events is too small to be 
measured, the second scenario (17) looks rather promis- 
ing. Our conceptually simple approach allows us to in- 
corporate the detector set-up into the final formulas for 
the partial cross sections relevant to experimental obser- 
vation at the LHC. Essentially, we conclude that in the 
first scenario (15), bound-free pair production takes place 
but most of the positrons escape the detector; they are 
not sufficiently deflected out of the beam line to be within 
the acceptance range of the LHC detectors. 

In the present work we have restricted our theoreti- 
cal analysis to a single e + e~ production. In the ultra- 
relativistic regime, however, the ion-ion collisions may 
result in the creation of a few electron-positron or even 
muon-antimuon pairs. Such a multiple-pair (and heavier 
lepton) production now attracts considerable attention as 
a tool for exploring the quantum electrodynamics in ex- 
tremely strong electromagnetic fields [34, 35]. In heavy ion 



colliders, these processes again can be explored by measur- 
ing the emitted positrons and residual few-electron ions in 
coincidence. Theoretical predictions, required for prepar- 
ing and analyzing such coincidence experiments, may be 
naturally obtained within the framework of the EPA- 
Sauter theory. Investigations along this line are currently 
underway and will be presented in a future paper. 
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